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We investigate different physical properties of a spin imbalanced fermionic superfluid described 
by an attractive Hubbard model in presence of a harmonic trap. To characterize the ground state 
of such a system, we compute various correlation functions, such as pairing, density and current 
correlations To underscore the effect of the trap, we compute these quantities as a function of 
strength of the trapping potential. A distinguishing feature of a superfluid state with population 
imbalance, is the appearance of a phase with spatially modulating gap function, called the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) phase which is present for both with and without trap. Trapping 
effects induce lowering of particle density and yields an empty system in the limit of large trapping 
depths. 



PACS numbers: 74.20.Fg 



I. INTRODUCTION 



The experimental realization of Bose-Einstein conden- 
sation in dilute alkali gases [IHH has triggered vast explo- 
ration of the field of ultracold, degenerate gas of fermionic 
atoms in presence of a harmonic trap. Achieving degen- 
eracy in a system of fermionic atoms is a well known 
hurdle as s-wave collisions between fermions in the same 
state are blocked by the exclusion principle (scattering 
in higher angular momentum channels, such as p-wave, 
are anyway weak). However researchers continued to ac- 
cess the degeneracy limit and the rare earths are recent 
additions on the list 0, @|- In most atomic gas experi- 
ments, the experimental challenge of producing ultracold 
temperatures is achieved by confining the atoms using a 
trapping potential. In fact, the velocity or momentum 
distribution of the gas is obtained via a technique called 
'time of flight' which involves a sudden turn off of the 
trapping potential and hence wait for an expansion time. 
The trapping potential is produced due to interaction 
between a gaussian laser beam and the induced atomic 
dipole moment. At the bottom of the trap, the gaussian 
profile can pretty much be approximated by a harmonic 
function, except for very shallow trap depths. 

In this paper, we shall incorporate the optical trap 
depth effects by a harmonic trap. So our system com- 
prises of a population imbalanced superfluid in the pres- 
ence of a harmonic trap which has a minimum (deep- 
est) at the center of the (two-dimensional) lattice and 
gradually increases towards the edges [a, |9(. Further we 
shall focus on the effect of the trap on properties of the 
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ground state by reviewing different correlation functions 
that characterize the state. These are sure to provide 
important insights on the nature of the polarized Fermi 
gas. 

A convenient technique to emphasize the effect of the 
trap is to examine the scenario in absence of a trap and 
compare it with those for various trap depths. Further, 
since the presence of a harmonic trap forces the system, 
which is a two-dimensional superfluid of Fermi atoms, 
to loose translational invariance, it is prudent to investi- 
gate correlations in the real space. In this paper, we in- 
clude several of them such as, gap parameter, off-diagonal 
long range order, pair-pair, density-density and current- 
current correlations to study the ground state of the pop- 
ulation imbalanced Fermi gas in presence of a harmonic 
trap. The above list may not be exhaustive, but are 
certainly useful to characterize a state that incorporates 
physics of superfluidity in presence of combined effects of 
population imbalance and harmonic trap. 

The results that we obtain have two special features. 
A phase with a finite momentum Cooper pairing and 
characterized by a spatially modulating order parame- 
ter, exists for moderate values of polarization of the par- 
ticipating species of fermions. This state, called as the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, persists 
for both with and without trap[10j. However large im- 
balance gives way to a fully polarized Fermi liquidfllj]. 
Another interesting artifact of the presence of the trap 
occurs in the form of loss in density of the constituent 
particles of the system. This is inevitable since the chem- 
ical potential is required to be kept fixed to attain a pop- 
ulation imbalance of the participating species and thus is 
crucial for observing the FFLO phase. Hence in the limit 
of large trap depth, we obtain an almost empty system. 
This almost empty system has some features that mimic 
an insulating phase, such as a vanishing superconducting 
gap, an increasing spectral gap and a dramatic reduction 
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in mobility of the particles as a function of the trapping 
potential. It should be noted that a gradual emergence 
of such an empty system with increasing trapping effects 
can not be termed as a transition to an insulating state 
which are otherwise present in the repulsive version of 
the Hubbard model[12;]. Such transition have been ex- 
perimentally realized in case of bosonic atoms [l3| that 
are supported by theoretical studies [14|. 

The study of population imbalanced gases of fermionic 
atoms in presence of a harmonic trap has a reasonable 
history [8l. l9l. 1 1 5Ml 7| . Density distribution of the fermionic 
atoms in the trap has been a central issue and a phase 
separation is predicted with the paired atoms confined at 
the center of the trap, while the unpaired ones find place 
towards the periphery [H However a concomitant 
tuning of the spin imbalance and trap depth was lack- 
ing. This paper fills the void and explores the interesting 
physics associated with it. 

We organize our paper in the following fashion. In 
order to fix notations, we provide a description of the 
attractive Hubbard model in presence of a harmonic 
trap and the standard mean field Bogoliubov de Genncs 
(BdG) theory which are used here. In the next section, 
we present results for different correlation functions in the 
presence of a trap and compare with the results when the 
trapping effects are switched off. For a range of popula- 
tion imbalance caused by a constant Zeeman field in our 
model, signature of FFLO phase is observed. The FFLO 
phase is characterized by weaker correlations but has a 
robust presence both with and without trap effects [1]. 
Further as a function of the trap depth, we demonstrate 
a gradual lowering of the particle density in table I. A 
probe into the nature of the almost empty state yields 
a zero superconducting gap, an increasing spectral gap 
and a reduced average kinetic energy of the particles as 
the strength of confinement is enhanced. These features 
may indicate emergence of an insulating-like state, how- 
ever when considered along with a gradual vanishing of 
particle density, it rules out the possibility of a transi- 
tion from a supcrfluid to an insulator. We conclude with 
a brief discussion of our results. 



II. MODEL AND FORMALISM 

In crystal lattices, where electrons are the main carri- 
ers that contribute to transport properties of materials, 
obtaining a model Hamiltonian that describes the phys- 
ical properties satisfactorily is a difficult task owing to a 
highly complicated band structure. However a gaseous 
system of fermionic atoms in a confining potential is a 
much neater realization of the simplest model that in- 
corporates strong correlations between atoms at short 
distances. An attractive version of the two-dimensional 
Hubbard model in which the atoms interact attractively 
via the on-site interaction, U in presence of a constant 
Zeeman field, h and a harmonic trapping potential at site 
Tii Vi is written as, 



«= "* E ( C U CT +H.c.) (1) 

- 1^1 E( n *t ~ \)( n n " \) + J2( Vi ~ v + <jh ) n ^ 

fi denotes the chemical potential, t is the hopping matrix 
element between the nearest neighbours of a two dimen- 
sional square lattice. The creation (annihilation) oper- 
ator for fermionic atoms corresponding to spin a is c\ a 
(cia). The excess of one species of atoms (say with spin- 
f) over another is controlled by the magnetic field h (or 
equivalently an effective chemical potential, jj = ah — /i). 
Vi is assumed to be of the form, 

Vi = Vo(n - r ) 2 (2) 

where Vq is the strength of the trapping potential and ro 
is the position where the center of the trap lies and is 
located at the center of the lattice. Thus the potential 
is minimum (deepest) at the center of the lattice and is 
maximum (shallow) at the edges. All of U, h, /j, and Vb 
are expressed in units of hopping strength, t which is typ- 
ically of the order of an eV. A mean field decoupling of 
the interaction term in Eq. (fTJ) yields the effective Hamil- 
tonian of the form, 

ij,a i 

here A$ = — |Z7|(cjj,Cj-|-} is the gap parameter for the 
fermionic superfluid. Hi jo — ~t^i±ij + {Vi— fi — U Sntg- + 
ah)Sij where 8n lS = n lS - 1/2 with (n la ) = (c\ a c ia ) and 
a = —a. 

Eq. ([3]) is hence diagonalized using Bogoliubov trans- 
formation which yields, 

f % A, \ fu n (r,.n = Ju n { Ti ) \ 

\ a, -uy) \ v n {*i) ) n V w «( r *)/ 

where u n (rj) and v n {vi) are the BdG eigenvectors satis- 
fying ^„[iij(r,) + v*(ri)] = 1 for all r, and E n are the 
eigenvalues. 

The gap parameter and density (and hence magneti- 
zation, TOj = (riif) — (nil), see Eq. ©) in terms of the 
eigenvectors u n (vi) and v n (vi) at a temperature T are 
given by, 

A, - - \u\J2K(*i)<(n)f(Ent) (5) 

n 

- u n (ri)v*(ri)f(-E ni ,)] 

M = [\ U n(*i)\ 2 f(Ena) + |t>„(r i )| 2 /(--Bnff)] (6) 
n 

where f(E na ) is the Fermi distribution function. A, 
and (niu) are obtained self-consistently from Eq. J5| and 
Eq. ([5]) at each lattice site. 
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It may be noted that a number of self-consistent solu- 
tions may exist for the gap parameter, A, corresponding 
to one set of parameters with different initial guesses. 
The winner among these will be decided by computing 
the free energies, T computed with respect to the free en- 
ergy of vacuum (at zero temperature) and is given by[20lj . 



J 7 — E E n c 



f(E na )-J2\v n (r t )\ 2 



(7) 
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We further compute higher order correlations which 
are likely to be more illustrative in characterizing the 
phase. We begin with off-diagonal long range order 
(ODLRO) which characterizes the superfluid phase and 
can be useful in bringing out distinction between a zero 
net center of mass momentum pairing (BCS) and a fi- 
nite momentum pairing states (FFLO). The ODLRO or- 
der parameter, A op is defined by the long distance be- 



haviour of the correlation 



for |r ? - — r,-| — > oo 



-»• A^ P /|C/| 2 
|2l|. The other correlation functions 



which are of physical significance are the pair-pair (CV, — 

( c lt c !i c Jt c jV)' density-density (K?f = {c\ a c ia J ja ,c J<7 >)) 
and paramagnetic current-current correlation functions, 
Ay 3 = (j a {ri)jfs(rj)) where j a is the component of the 
paramagnetic current density operator at site Ti given 



by j a (n) = itJ2 c 



Here r, and 



Tj denote lattice sites 



i+a,a C i,& ~ C i,a C i+a,<T 

In our numerical computation, 
a pair (say for Cij) is located at a fixed site and the 
probability of another pair at site rj, where rj can be 
any other site of the lattice, is calculated. Similar nota- 
tion is followed for K°° and A^ 3 as well. It is useful 
to mention here that while we have considered site rj to 
be along any arbitrary direction with respect to site rj, 
however results (shown in Fig. 6-8) are presented for rj 
to be along the length (x-axis) of the lattice. The rea- 
son is the following. While the direction along which rj 
is chosen is immaterial for the homogeneous BCS phase, 
for the FFLO phase, the preferred direction (for the case 
of without trap) is set by the propagation of the modu- 
lating gap parameter, which is along the length (x-axis) 
of the lattice (Fig. 2(c)). To retain uniformity, we follow 
the same for the trapped case, where modulation exists 
along a ring away from the trap center. 

Finally, we comment on the choice of parameters. 
We perform our studies for an inter-particle attraction 
strength \U\ = 2.5t which may be considered to be weak, 
and hence suitable as we propose a BCS superfluid as 
a starting point. It should be noted that further lower 
values of \U\ yield a vanishing gap parameter for any 
reasonable value of density. We have chosen \i to be 
— 0.5i. All quantities are calculated at zero tempera- 
ture, where the Fermi function has been taken as unity. 
We have also considered other values for density, rang- 
ing from moderate to large (very small densities do not 



yield a stable superfluid phase) and different (weak) in- 
teraction strengths, U, however we have not noted any 
qualitative change in behaviour of the correlation func- 
tions which are on focus in this paper. Further, it may 
be noted that the gap parameter undergoes a one di- 
mensional modulation[22] with a period that is commen- 
surate with the lattice size in the absence of confining 
potential. Thus a two dimensional lattice of size 32 x 16 
has been considered so that more periods of modulation 
of the gap parameter can be accommodated. We have 
also checked that the results for a 32 x 32 lattice are unal- 
tered with regard to the modulation of the gap parameter 
and qualitative behaviour of the correlation functions. In 
case of presence of the external harmonic potential, the 
dimension of the lattice is chosen to be 24 x 24 with open 
boundary conditions and the trap center is located at the 
center of the lattice. 



III. RESULTS 

To characterize the spin polarized phase in presence of 
a harmonic trap, we have computed various correlation 
functions in the real space for various values of popula- 
tion imbalance and trap depth. The different correlations 
that we consider in this work are gap parameter, mag- 
netization, off-diagonal long range order (ODLRO), pair, 
density and current correlations. The list is not exhaus- 
tive but will serve our purpose in providing insights on 
the nature of the population imbalanced superfluid in the 
superposed environment of a harmonic trap. 

Before we proceed with the discussion on the results 
for various correlation functions, we wish to distinguish 
between the three phases, viz. the homogeneous BCS 
phase for low values of magnetic field (i.e. low popula- 
tion imbalance), a spatially inhomogeneous FFLO phase 
for moderate field strengths and finally a normal phase at 
higher magnetic fields. Thus to fix our input parameters, 
we present a schematic diagram of the BCS and FFLO 
phases, both in the absence and presence of trap in Fig. 1. 
The FFLO phase is marked by a lower and a upper crit- 
ical magnetic field values, h Cl and h C2 respectively and 
lie intermediate to the BCS and normal phases. These 
boundaries are obtained from the behaviour of the gap 
parameter obtained solving Eq. (J5J) . Fig. 1 shows that the 
values for h Cl and h C2 are 0.35i and 0.5t in the absence 
of trap whereas it shifts to h Cl = 0.25t and h C2 = 0.4t 
when the trapping potential is switched on for \U\ — 2.5t 
and \i = — 0.5t. Hence all our results are presented for 
two values of magnetic field, h which denote two differ- 
ent values of population imbalance. One being in the 
homogeneous BCS phase (h < h Cl ) and another in the 
FFLO phase (h Cl < h < h C2 ). From Fig. 1, we choose 
h = O.lt and h = 0.5t corresponding to no trap and 
h = and h = 0.3i in presence of the harmonic trap 
for BCS and FFLO phases respectively. Moreover, we 
have restricted the demonstration and discussion of the 
results to one value of strength of the trapping potential 
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viz. Vq = 0.0I6i. The representative value considered 
here, approximately denotes a point where there possi- 
bly exists a very shallow minimum in the spectral gap. 
We have investigated the correlations for other values of 
trapping depths in both the superfiuid (small Vo) and 
strongly confined regimes (Vo ~ 0.2i, see Figs. 5 and 9 
and related discussion). 



BCS 



FF L O 



Normal 



FIG. 1: A schematic representation of the FFLO phase is 
shown as a function of the magnetic field both in the absence 
and presence of trap for \U\ = 2.5t and fi — —0.5t (same for 
all plots) . This phase is intermediate to a BCS phase (homo- 
geneous Aj) and a normal phase (Aj = 0). The boundaries of 
the FFLO phase are marked by h ci and h C2 - The solid lines 
are the boundaries in the absence of trap whereas the broken 
lines mark the boundary in the presence of the trapping po- 
tential of strength, Vo = 0.016t. Same numerical value of Vo 
is chosen for Figs. 2-4 and Figs. 6-8. Also all the quantities 
are in units of hopping integral, t. 

We now present our results for the gap parameter, A, 
in presence of a harmonic confinement and compare with 
the corresponding behaviour in absence of the trap. It 
may be noted that Aj is homogeneous for low values of 
magnetic fields (population imbalance being small) in the 
absence of trap, whereas the trapping effects lead to a 
maximum of Aj at the trap center, hence decreases and 
finally vanishes away from the center (see Fig. 2 (a) 
and (b)). Such contrasting behaviour is caused by the 
trapping effects of the harmonic potential that results in 
the accumulation of atomic densities around the center 
thereby leading to such an inhomogeneous distribution 
of Aj in real space. Next we compare Aj for the modu- 
lating (FFLO) phase. Aj modulates along one direction 
(x-axis) of the lattice without the trap while it modu- 
lates along the radial direction in the presence of trap, as 
shown in Fig. 2 (c) and (d). That is, far away from the 
trap center, the gap parameter undergoes a sign change 
confirming the presence of the FFLO phase. 

We proceed to discuss results on the local magnetiza- 
tion, TOj (= (riif) — {riii)) presented in Fig. 3. While m; 
for the BCS phase, is zero understandably (being non- 
magnetic due to lack of unpaired particles) across the 
lattice (see Fig. 3(a) and (b)), it modulates with the pe- 
riod half of that of Aj in the FFLO phase when trapping 
effects are not included (Fig. 3(c)). The large values of 
magnetization occurs at nodal lines with broken pairs, 



whereas lattice sites with finite Aj correspond to weak 
magnetization, thereby leading to a phase difference be- 
tween these two quantities. The magnetization profile 
in the presence of trap, exhibits a bimodal structure (as 
shown in Fig. 3(d)) since the magnetization is non-zero 
around the ring like nodal line. In other words, a small 
number of particles are squeezed into the inner core of 
the harmonic trap due to pair formation in the super- 
fluid phase, while the unpaired (majority) carriers are 
pushed out to outside of the core. 

We now plot the ODLRO order parameter, A op as a 
function of magnetic field without the harmonic trapping 
in Fig. 4. It shows sharp drop of A op at the onset of the 
FFLO phase at h Cl — 0.35i, thereby indicating weak- 
ened superconducting correlations in the FFLO phase 
(see Fig. 1). A op is also computed in the presence of 
trap but it is seen at two fixed values of magnetic field 
viz. h = and h — 0.M for various values of trapping 
potential, Vo (shown in Fig. 5). Our results clearly in- 
dicate a higher value of A op for h = (no population 
imbalance) than that for h = 0.3t, but is suggestive of a 
rapid decay for both with increasing confinement of the 
carriers. 

Next, we compute pair-pair (Cy), density-density 

(Kfj 7 ) and current-current (Ay ) correlation functions, 
the utilit y o f which in experiments are illustrated in 
literature |23l - [25| . All these quantities point towards a 
spatially homogeneous (translationally invariant) ground 
state for the BCS phase, while the translational invari- 
ance is broken in the FFLO phase resulting in spatial 
inhomogeneities of these correlation functions. The cor- 
relation functions mentioned above need individual at- 
tention. 

We begin with the pair-pair correlation function, Cy 
(defined in the previous section), the real space scan of 
which taken along the length of the lattice (along the di- 
rection in which Aj modulates), from one end to another 
is presented in Fig. 6(a) and (b). Cy in the BCS phase 
(h = 0. It) and without an external confinement, is higher 
(approximately ten times) than the value in the FFLO 
phase (h = 0.5i). This supports weakening of the super- 
conducting correlations with increasing magnetic field as 
also is evident from Figs. 2(a) and (c). The behaviour of 
Cy in the presence of trap are interesting in the following 
sense. It shows a maximum value at trap center and falls 
off in a near symmetric fashion away from the trap cen- 
ter. It also indicates towards higher correlations in BCS 
(h = 0) which are considerably reduced (approximately 
to one-third) in the FFLO phase (h = 0.3t). 

The reason behind Cy showing a hump-like behaviour 
is as follows. The probability of finding a pair in the 
center of the trap is higher as the underlying trapping 
potential is minimum (deep) at the center. However, it 
drops at the edges where the trapping potential is high 
(shallow). Additionally, Cy drops significantly as the 
external trapping potential is switched on, suggesting the 
fact that the trapping potential assists an imbalanced 
superfiuid phase to have weakened pair-pair correlations. 
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FIG. 2: (Colour online) Local pairing amplitude, Ai in BCS ((a) and (b)) and modulating FFLO ((c) and (d)) phases are 
shown. Figs, (a) and (c) are for cases without trap while (b) and (d) are in presence of trap. The system size is 32 x 16 in the 
absence of trap and 24 x 24 when the trapping potential is present and are same for all our results. 



Hence we discuss density-density correlation function, 
K%? . In absence of the trap, there is no significant dif- 
ference between the behaviour of Kff in the BCS and 
FFLO phases (see Fig. 7(a)) with both remaining con- 
stant throughout the lattice. However, the correlations 
show a hump-like feature as the trapping potential is 
turned on (Fig. 7(b)). Though the shape of the plot is 
similar to that of Cij, the qualitative behaviour in both 

the phases is reversed i.e. Kf° for the FFLO phase 
(h = Q.3t) is higher than the BCS phase (h — 0). It may 
be noted that we have considered all a, a =t, i- However 
the results are only shown for a = a = J [. Other choices 
yield qualitatively same results. A possible explanation 
for the rise in correlations between different species (t 
and l) of atoms for the FFLO phase can be explained in 
terms of the increase in number of unpaired particles with 
the increase in the magnetic field. Thus the correlations 



are higher in the FFLO phase (with more number of un- 
paired carriers) as compared to the BCS phase where all 

carriers are paired. Further, a comparison of Kf? in the 
presence of trap to the one without trap shows remark- 
able reduction in the amplitude in the presence of trap. 
These correlations are particularly useful in time of flight 
experiments in which the trapping potential is suddenly 
switched off, resulting in free expansion of the atomic gas 
and thus enable realization of the velocity or momentum 
profile [24]. These techniques have been used to detect 
the nature of the condensate e.g. the superfluidity of the 
atoms constituting the condensate etc|26j|. 

Finally, we analyze the behaviour of the paramag- 
netic current-current correlations (with a = j3 = x), i.e. 
A^j in an environment of harmonic confinement. Af? 
rapidly decays to zero in the BCS phase, while it mod- 
ulates across the lattice and remains finite in the FFLO 
phase, irrespective of the value of the trapping potential 




h = 0.3f 



FIG. 3: (Colour online) Local magnetization, m% is shown in BCS (a) and (b) and FFLO (c) and (d) phases both in the absence 
and presence of the confining potential. 





0.01 0.015 
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FIG. 4: The ODLRO order parameter, Aop (in units of t) is 
shown as a function of magnetic field, h. The dotted line is a 
guide to the eye and used in some of the subsequent plots. 



FIG. 5: The ODLRO order parameter, Aop (in units of t) is 
shown as a function of trapping potential, Vo for h — (bold 
line) and h — 0.3t (points and dashed line). 



as shown in Fig. 8(a) and (b). The modulation in real possibly because of the involvement of more than one 
space does not have a uniform profile as Aj (Fig. 2(c)), momentum vector in Af? Nevertheless, we still have the 
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FIG. 6: Pair-pair correlation function, Cy is shown across the lattice in BCS and FFLO phases without trap (a) and compared 
with the case when the trapping effects are invoked (b). 
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FIG. 7: Density-density correlation function for a = a ="f i.e. KjJ is shown across the lattice in BCS and FFLO phases 
without trap (a) and with trap (b). 



response of the system to an external perturbation, to be 
more in the FFLO phase than that in the BCS phase and 
thus in turn reiterates the presence of stronger supercon- 
ducting correlations in the homogeneous BCS phase. We 
have recorded higher values for Af? for the case where 
no trap effects are included than those in presence of the 
trap. The likely reason must be the reduced mobility 
of the charge carriers, thereby causing weakened current 
vectors in presence of the trap due to localization effects. 

In the following discussion, we describe the scenario 
corresponding to strong confining effects. As the trap 
is made deeper, in order to keep the chemical potential 
fixed (needed to obtain a population imbalanced state 
and FFLO phase to be realized), the particle density is 
lowered. So the system looses particles and heads to- 
wards an empty state with only a few particles remaining 
at the core of the potential as trapping effects are made 
stronger (table I). 



In order to obtain a detailed understanding of what 
happens at larger trap depths, we compute certain phys- 
ical properties that characterize the state. The gap pa- 
rameter (Aj) vanishes for this scantily populated system 
(not shown here). The spectral gap, E gap which is de- 
fined as the difference in energy between the ground state 
and lowest lying excited state of the BdG spectrum as ob- 
tained by numerically solving Eq. @, shows an increase 
with increasing Vo- There is possibly a very shallow min- 
imum that exists around Vo ~ 0.016t (Fig. 9(a)) as seen 
in our numerical computation, however this feature is not 
clear from Fig. 9(a). Since an increasing E gap suggests 
the onset of an insulating behaviour (27|, we explored the 
average kinetic energy of the particles (Fig. 9(b)) as a 
measure of their mobilities and find a rapid decay as Vb 
is increased. 

All these results put together may suggest of the onset 
of an insulating phase at large trap depths, however when 
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FIG. 8: Current-current correlation for a = f3 = x i.e. Af/ is shown across the lattice in BCS and FFLO phases without trap 
(a) and with trap (b). 




0.12 0.15 




1 0.15 y 0.2 



FIG. 9: (a) The spectral gap E gap is shown as a function of increasing trapping potential Vo. (b) The average kinetic energy 
of the charge carriers, (KE) is plotted as a function of Vo- Both plots are presented at h = OAt. 



taken together with the fact that particles are gradually 
leaking out (see table I), they appear to be artifacts of 
emptying of the system. 

IV. CONCLUSIONS 

In summary, we have explored the ground state of 
a population imbalanccd Fermi gas on a two dimen- 
sional lattice in presence of a harmonic confinement. 
We adopted a weak coupling s-wave superconductor de- 
scribed by an attractive Hubbard model in presence of 
imbalanced spin species induced by a Zeeman field. The 
ground state is inhomogeneous for a range of magnetic 
field strengths which is characterized by a modulating or- 
der parameter, i.e. the FFLO phase. The exotic FFLO 
phase is obtained both in presence and absence of the 
trap. 



To aid our understanding of the trapping effects on the 
system, we computed various correlation functions e.g. 
off-diagonal long range order, pair-pair, density-density 
and current-current correlations in real space and pro- 
vided adequate comparisons with the scenario when the 
confining potential is absent. The results are indicative 
of weakened correlations as population imbalance is in- 
creased and confining effects are invoked. 



Larger confining effects induce loss of particles from 
the system (chemical potential being fixed to its original 
value). At large values of Vo, the system becomes almost 
empty of particles. The physical properties of this empty 
system behaves similar to that of an insulator. However 
such signatures should not be considered indicative of 
a transition from a superfluid to an insulating phase as 
they are artifacts of a near empty system. 



9 



TABLE I: The total density of particles, (n)(= ^Y, CT ( n i<r)) is 
shown as a function of increasing trapping strength, Vo (in 
units of t). The density drops significantly as the trapping 
strength is made stronger leading to an empty system at large 
values of Vo- 
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0.00 


0.6361 


0.016 


0.2135 


0.05 
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